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Abstract: 

We discuss the concepts and methodology to implement an experiment probing directly 

f^ '_ Time Reversal (T) non-invariance, without any experimental connection to CP violation, 

by the exchange of in and out states. The idea relies on the B°B° entanglement and 

decay time information available at B factories. The flavor or CP tag of the state of 
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the still living neutral meson by the first decay of its orthogonal partner overcomes the 



J> . problem of irreversibility for unstable systems, which prevents direct tests of T with in- 



coherent particle states. T violation in the time evolution between the two decays means 
experimentally a difference between the intensities for the time-ordered (£ + X, J/i/jK s ) and 
(J/ipK L , £~X) decays, and three other independent asymmetries. The proposed strategy 
has been applied to simulated data samples of similar size and features to those currently 
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CNJ ■ available, from which we estimate the significance of the expected discovery to reach many 

standard deviations. 
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1 Introduction 

Violation of the CP symmetry has been observed in K and B meson decays [1—4]. In the 
Standard Model (SM) the mechanism of CP violation in weak interactions arises from the 
joint effect of the three non- vanishing mixing angles and the single irreducible phase in 
the 3-family Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix. The existence of this 
matrix conveys the fact that the quarks that participate in weak processes are a linear 
combination of mass eigenstates. This mechanism has been confirmed to be the dominant 
by experiments probing CP violation, particularly with studies involving B decays [5]. 

In the context of local quantum field theories with Lorentz invariance and Hermiticity, 
the CPT theorem [6] ensures a theoretical constraint between the CP and T symmetries. 
Although all present tests of CPT invariance confirm the validity of this symmetry, par- 
ticularly in the neutral- kaon system where there are strong limits [7, 8], the theoretical 
connection between CP and T does not imply an experimental identity, except for pro- 
cesses which are identical under CPT transformation. Moreover, it is worthwhile to search 
for direct evidence of T non-invar iance, with neither experimental nor theoretical connec- 
tion to CP violation and CPT invariance. There is at present no existing result that clearly 
shows T violation in this sense [9]. 

In the case of transition processes, due to the antiunitarity of the operator implement- 
ing the symmetry, T invariance requires that the likelihood for reaction in — > out equals 
that of out — > in once the initial configurations, namely in in one case and out in the other, 
have been precisely prepared. However, the likelihood of the time reversed version of the 
process to happen is very low or impractical in some cases, like for unstable systems. This 
explains why T is much more difficult to study directly than P, C, and CP. 

Since the SM is CPT invariant, it predicts T- violating effects in parallel to each CP- 
violation effect that arises due to the interference of amplitudes with different weak phases. 
These may appear in three different ways: T violation in decay; T violation in the mixing 
of neutral states; and T violation that arises from the interference between decay with and 
without mixing. 

T violation matched to CP violation in decay has not been observed, because the dif- 
ficulties of the preparation of the time reversed decay process. Let us see the example [10] 
of the rare weak decay of neutral- B mesons to K + ir~, for which direct CP violation is 
well established [4] (different decay rates R\ and R2 for B° — > K + ir~ and B° — > K~ir + , 
respectively) due to the ability of B factories to produce hundreds of millions of B mesons. 
T violation as implied by this result, combined with CPT invariance, tells us that the rates 
for the inverse processes K + n~ — > B° and K~tt + — > B° should be R2 and R±, respectively. 
However, there is little chance to measure these inverse rates and check directly this predic- 
tion since the weak interaction production mechanism is highly suppressed (B° — > K + ir~ 
branching ratio of order 10 -5 ) and the strong interaction would completely swamp the 
feeble weak process. 

T violation associated to CP violation in the mixing has been experimentally analyzed 
in K [11] and B mesons [12]. Here one looks whether the rate for a neutral- K (B) meson 
tagged at its production as K° (B°) and identified afterwards as K° (B°) is equal to the 



rate for the neutral particle tagged at its production as K° (B°) and identified later as 
K° (B°). Any difference in this case is both CP and T violating, because CP and T 
are experimentally identical for this process. The experimental results for kaons yielded a 
T-violating difference in these rates. Such a difference is proportional not only to the T- 
violating term of the K°K° matrix that defines the mass eigenstates in terms of the flavor 
eigenstates, but also to the width difference Ar between the two mass eigenstates, thus 
T non-invariance would not be present in the limit Ar — > 0. Therefore, this asymmetry 
shows T violation proportional to Ar, time independent, experimentally identical to CP 
violation, thus it is not an independent T non-invariance test as one might like [9]. In the 
B°B° system no asymmetry has been yet found, as expected within the SM since in this 
case Ar almost vanishes [13]. 

The largest CP-violating asymmetry in Nature has been found between the rate for 
B° — > J/ipK s (and other similar CP-odd ccK s final states, e.g., ip(2S)K s , XciK s , or 
the CP-even J/t/)K L final state) and the CP-conjugate rate for B° to decay to the same 
CP-eigenstate, which is originated in the interference between the time-dependent decay 
amplitudes with and without mixing [3, 5]. Here, CP violation arises because the mixing 
phase, i.e., the relative phase between the complex parameters defining the mass eigenstates 
in terms of the flavor eigenstates [14] minus the relative phase of the ratio of the amplitude 
for the decay and its CP conjugate, does not vanish. In the SM this phase difference is 2/3, 
where /3 is the angle between the V c dV*b an< ^ ^tdYtb sides of the CKM unitarity triangle. 
At B factories the experimental studies are performed using the two B mesons produced 
in the antisymmetric coherent state from the T(4S) decay. One neutral- P> meson decays 
into a definite flavor state and the other is reconstructed in the CP-eigenstate final state 
of interest with a given decay time difference At = tcp — ^flavor- The measured asymmetry 
is large, proportional to sin 2/3 « 0.7, a time-odd dependent function that reverses sign 
between B° and B° tagged events, and between B° — > J/ipK s and B° — > J/ipK L events 
with the same flavor tag. This is the observed CP-violating effect. However, since there is 
no reversal of in and out states, this time asymmetry cannot be interpreted as genuine T 
violation. 

In this paper we propose and describe a methodology that makes use of the Einstein- 
Podolsky-Rosen (EPR) entanglement [15] available at current B factories to overcome the 
problem of irreversibility and perform a direct observation of Time Reversal Violation. The 
method relies precisely on the possibility for preparing the quantum mechanical individual 
state of the neutral- P meson by the observation of particular decay channels of its orthogo- 
nal entangled partner, and studying the time evolution of the filtered state of the still living 
meson. This strategy allows the interchange in <->■ out states for a given process. Whereas 
the basic ideas have been presented previously [16] and scrutinized later [9, 10, 17], here 
we discuss for the first time the steps to implement these concepts into an actual experi- 
ment able to produce the desired results. In addition, using a simulation of data samples 
of similar size and properties to those currently available at B factories, we illustrate the 
application of the methodology and estimate the significance of the expected discovery. 
The proposed analysis is independent of the underlying theory, CPT invariant or not, for 
describing the relevant transitions. 



2 Time Reversal from entangled B°B° decays 

Bose statistics, for a B°B° system with C = — , requires these two mesons to be in an 
antisymmetric state of the system for any pair of orthogonal individual states. These can 
be either the B states projected by flavor, 

\i) = ±=lB\ tl )B\t 2 ) - 5°(ti)S°(fc)], (2.1) 

or the ones projected by CP final states [19], 

\i) = i=[ J B + (t 1 )B_(t 2 ) - B_(t 1 )S+(t 2 )], (2.2) 

where t\ and t% are the labels to specify the states "1" and "2" of each neutral- B meson 
by means of the time of its future decay, with At = ti — t\ > 0. The antisymmetry 
remains invariant with the time evolution, including mixing, before the first decay at t\. In 
Eq. (2.2), 2?_ is the neutral- .B state filtered by its decay to J/ipK + , K + being the neutral- 
K state filtered by its decay to 7T7T, and B + is orthogonal to -B_, not connected to J/ipK + . 
Similarly, in Eq. (2.1), B° is the neutral- .B state filtered by its decay to £~X, for example 
a semileptonic decay with a negatively charged lepton or a hadronic final state containing 
a D + or D* + meson. We note that treating J/ipK s , K s — > irir final states as J/ipK +) i.e., 
neglecting CP violation in the neutral-kaon system, introduces effects that are small and 
thus can be neglected. Similarly, treating hadronic final states containing DW+ mesons as 
B°, i.e., neglecting CP violation due to interference from doubly CKM-suppressed decays, 
introduces effects small enough that can be considered as corrections. The same applies to 
B + as the state connected to J/ifiK- and treated as J/ipK L , and B° as the state connected 
to hadronic final states containing D^*>~ mesons. The identification of B-, B + as those 
states filtered by the decay to CP eigenstates is illustrated in Appendix A. 

Therefore, in addition to the flavor tagging used in standard CP studies at B factory 
experiments [3, 18] we can apply a CP tagging [19] to one of the B mesons decaying into 
the CP-odd final state J/^K S , preparing the orthogonal B meson as B + state at the initial 
time t\. Afterwards it decays at ti and is reconstructed in the flavor (B° or B°) final state 
of interest (£ + X or £~X, respectively). This combination of flavor and CP tags allows to 
filter initial and final states to compare, for example, the rate for a B° evolving to B— 
(£ + X decay product of the B partner first, J/i(jK s final state later) with the rate for a 
B^ evolving to B° (J/ipK L decay product of the B partner first, £~X final state later). 
The relation between these T-conjugated transitions and the reconstructed final states in 
the experimental B factory scheme is illustrated in Fig. 1. The comparison between the 
time evolution of the neutral- B meson from its preparation as B° until its decay as B_, 
B°(t\) — > B-fa), and its T transformed B-(t\) — > B (t2) is our proposed test of Time 
Reversal symmetry. 

There are other three independent comparisons between T-conjugated processes, as 
summarized in Table 1. A non- vanishing asymmetry in the rates for any pair of T- 
conjugated transitions thus constitutes a direct and independent observation of T violation, 
in the sense discussed previously. 
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Figure 1. Sketch of two T-conjugate transitions in an experimental B factory scheme. The 
observation of the final states associated to the T-transformed transitions is divided into three well 
defined steps. We first observe the decay of one of the entangled B particles, produced in the T (4,5) 
decay, into a definite flavor (or a definite CP) decay products at fa , preparing the state of the other 
entangled B particle, which has not yet decayed at fa, into its orthogonal state. This tagged B 
meson state evolves in time to finally decay at <2 > fa into a CP (or a flavor) final state. It should be 
noted that T asymmetry is clearly different from the At (fa <H> £2) exchange and CP asymmetries. 
In fact, in the former we require to compare the reference transition B° — > _B_, flavor-tagged by 
e+X and decayed to J/^K S (1 + X,J/^jK s ), to the transition B_ -^ B°, CP-tagged by J/ipK L and 
decayed to £~X (J/ipK L ,£~X), whereas for At and CP asymmetries the reference decay products 
must be compared to (J/ipK s ,£ + X) and (£~X,J/ipK s ), respectively. 



We can also apply this methodology for similar tests of CP violation and CPT in- 
variance, providing a proof that T non-invariance is compensated by CP violation. Ta- 
bles 2 and 3 summarize all the possible comparisons of CP- and CPT-conjugated transi- 
tions, together with the corresponding final states. As anticipated, the transitions involved 
in the experimental tests of CP and T symmetries are different. For example, a test of 
CP symmetry can be done with the J/ijjK s final state only. On the contrary, a test of 
T invariance necessarily involves both J/tfjK s and J/ipK L final states. Furthermore, one 
may check that none of all comparisons between T-, CP-, or CPT-conjugated transitions 
in Tables 1, 2, and 3, respectively, corresponds to exchange of t\ and t<i- 



Reference 


T-conjugate 


Transition Final state 


Transition 


Final state 


B° -> B^ (£+X, J/^K S ) 


£_ ^B° 


(j/1>K L ,e-x) 


B + ^B° (J/^K S ,£+X) 


B°^B+ 


{l-X,Jj^K h ) 


B°^B+ (£+X,J/^K L ) 


B + ^B° 


(j/^k s j-x) 


B_^B° (J/iPK L ,£ + X) 


B° -> B_ 


(£-X,J/i,K s ) 



Table 1. Possible comparisons between T-conjugated transitions and the associated time-ordered 
decay products in the experimental B factory scheme. 



Reference 


CP- 


conjugate 


Transition Final state 


Transition 


Final state 


B° -»• B_ (£+X, J/ipK s ) 


B° -»• £_ 


(£-X,J/i;K s ) 


B + ^B° (J/^K S ,£ + X) 


B + ^B° 


(J/1>K 3 ,£-X) 


B°^B+ (£+X,J/^K L ) 


B°^B+ 


(£~X,J/^K L ) 


B_^B° (J/ipK L ,£ + X) 


£_ ^B° 


(J/rPK L ,£~X) 



Table 2. Possible comparisons between CP-conjugated transitions and the associated time-ordered 
decay products in the experimental B factory scheme. 



Reference 


CPT- 


conjugate 


Transition Final state 


Transition 


Final state 


B° -»• B^ (£+X, J/rpK s ) 


B_ -»• B° 


(j/iPK L ,e+x) 


B + ^B° (J/^K S ,£+X) 


B° ^ B + 


(£+X,J/tPK l ) 


B° -> £_ (£~X, J/ipKg) 


B_^B° 


(J/4>K L ,£-X) 


B+^B° (J/tJjK s ,£-X) 


B° -> B + 


(£~X,J/rpK L ) 



Table 3. Possible comparisons between CPT-conjugated transitions and the associated time- 
ordered decay products in the experimental B factory scheme. 



3 Methodology and T-violating parameters 

We can now proceed to a partition of the complete set of final states with definite flavor 
and CP content into 8 pairs, defined by the first decaying B at t\ and preparing the tagging 
state of the still living meson, i.e., B°, B°, B_, B + , as a function of At = ti — t\ > 0. Each 
of these 8 processes has a time-dependent decay rate g a JAt), where indexes a S {£ + ,£~} 
and /3 G {K S ,K L } run over the final states with definite flavor (£ + X, £~X) and CP 
eigenstates (J/i/jK s , J/ij)K L ), respectively, and the upper index + or — indicates if the 
decay to the flavor final state a occurred before or after to the CP-eigenstate final state (3. 
Thus a -\ > — replacement corresponds to At exchange, which means experimentally the 



exchange of the two decay products at £1 and £2- 

From only quantum mechanics each decay rate can be written as function of At 



g% B {A T ) oc e~ rAr { C± g cos(AmAr) + S± g sin(AmAr) + 



D±g cosh(ArAr) + £± sinh(ArAr) }, (3.1) 



where T is the average decay width, Am and Ar are the mass and width differences 
between the mass eigenstates, and C a g , S a g, D a g and E a g are generic coefficients. This 
construction makes no assumptions about neither CPT invariance nor CP or T violation. 
Assuming Ar = in the time dependence and renormalizing to the coefficient of the 
cosh(ArAr) term, Eq. (3.1) simplifies to 

g% 8 {Ar) oc e- rAr {l + C^cos(AmAr) + S^siii(AmAr)}. (3.2) 

The sine term in Eq. (3.2) results from the interference between amplitudes with and with- 
out mixing, whereas the cosine term arises from the interference between decay amplitudes 
with different weak and strong phases. Effects due to any small lifetime difference in the 
time dependence, and the renormalization of all coefficients to the cosh(ArAr) term for 
each subsample separately, introduce small corrections. 

It then follows that asymmetries in decay rates for any pair of T-conjugated transitions 
(Table 1) would be apparent through differences between their respective best fit S a g and 
C a o coefficients. For example, a significant difference between the S~[ + K and Sj_ K 
coefficients would imply observation of T violation. In general, a net difference of S a B or 
C aS parameters between two T-, CP-, or CPT-transformed processes (Tables 1, 2, or 3) 
would be a proof of T, CP, or CPT violation, respectively. 

The standard CP violation studies performed by the B factory experiments [18] extract 
a single set of best fit S and C coefficients, reversing the sign of S under At <->• —At, or 
B + f* B-, or B° f-> B° exchanges, and reversing the sign of C only under B° -H- B° 
exchange, 

(3.3) 



s 


_ 2/m(A) __ + + _ 

" 1 + |A| e+ > K s e-,K s °£+,k s °£-,k s 




q+ _ 0+ _ q— _ q— 
°£+,K L e-,K L £+,K L £~,K L > 


C 


_ 1 — A _ . _ + _ _ 
" 1 + |A| e+ > K s t-,K s ^1+,K S ^£~,K S 




e+,K L £-,k l £+,K L ^£-,k l - 



(3.4) 



This construction is valid under the assumptions of CPT invariance and Ar = [16]. If 
Ar = 0, automatically there is no CP violation in the mixing. Under these assumptions 
the At -H- —At (or equivalently t\ <R- £2) exchange, which is not a T-symmetry operation, 
becomes related to T (exchange of the in and out neutral- -B states) , and to CP (exchange 
of B° and B° states), CP f* T -H- At [20]. In other words, the resulting statement that 



particle <B- antiparticle invariance test is related to the exchange At -H- —At. In the SM 
A = e t2/3 , therefore 

S = sin2/3 (3.5) 

C = 0. 

As most of the currently available data on neutral- B mesons are well described by the SM, 
whether our proposed T-violating observables are also well accounted for in the SM or not, 
the direct observation of T violation would provide a proof of principle. 

It is now convenient to introduce the asymmetry parameters AS T , AC T (and similarly 
AS* CP , AC CP , and AS'pprp, and AC CPT ), defined in Table 4, as the difference of S a a, C a a 
coefficients between the corresponding symmetry-transformed processes for two reference 
At-exchanged transitions, for example Si K and S7 + K . Using these parameters rather 
than the S a g, C a o coefficients has the main advantage that the symmetry violation will 
be apparent through a non- vanishing value of any of the associated four parameters. In 
other words, if AS^ ^ or AS^ ^ or AC^ ^ or AC-p" / 0, then there is T violation, 
and similarly for CP and CPT symmetries. Table 5 summarizes the values and relations 
among the asymmetry parameters in the case of invariance under the three space-time 
discrete symmetry transformations. We would like to emphasize that the definition of 
these parameters requires to choose as reference two At-exchanged samples, in our case 
S £+ K since J/ij)K s events are cleaner and more efficiently reconstructed than J/ij)K L [18]. 

4 Simulation study 

We show in this section how the proposed methodology can be applied to an actual ex- 
periment to produce the desired results. The experimental samples are generated using a 
Monte Carlo simulation technique and are similar in size and properties to those currently 
available at B factories and used for standard CP- violation studies [18]. Thereby this study 
will provide an up to date estimate of the significance of the expected observation. 

4.1 Generation of samples 

Samples of events with a B° or B° state (i.e., £ + X or £~X) and a CP-eigenstate J/ipK s , 
ip(2S)K s , XciK s , denoted generally as ccK s , or J/ipK L , are generated taking into account 
their relative branching ratios, reconstruction efficiencies and misidentification rates. We 
consider about 8000 ccK s and 6000 J/ipK L events, with purities around 90% and 60%, 
respectively [18]. The signal component uses the probability density function (PDF) given 
in Eq. (3.2), with values for Am = 0.507 ps _1 , 1/r = 1.519 ps [8], and the coefficients 
S a o and C a a as given in Table 4 assuming the SM. The misidentification of CP final 
states is accounted for through a signal probability that depends on the two usual kine- 
matic variables at B factory experiments, the energy of the J/ij)K L and the mass of the 
ccK s candidates, including the different individual background sources [18]. Mistakes in 
the flavor identification (mistags) are included by modifying the time-dependent PDF as 
follows, 

h aA Ar ) °c (1 - Wa)^(Ar) + u a g±p(Ar), (4.1) 



Coefficient 


Assumed value 


Fit value 


AS+ = 


- <7~ <?+ 

£~,K L £+,K s 


-1.4 


-1.57 ±0.15 


A5t = 


- q+ o- 

J e-,K L °e+,K s 


1.4 


1.25 ±0.19 


AC+ = 


£-,K L t+,K s 


0.0 


-0.07 ±0.14 


AC T = 


£-,K L £+,K s 


0.0 


-0.09 ±0.14 


A5+p 


- q+ q+ 

£~,k s £+,k s 


-1.4 


-1.65 ±0.11 


AS" p 


l-,K s £+,K s 


1.4 


1.54 ±0.13 


AC+p 




0.0 


0.03 ±0.10 


ACcp 


£-,K s l+,K s 


0.0 


-0.09 ±0.10 


A5 CPT 


= S £+,K L ~ S £+,K S 


0.0 


-0.25 ±0.22 


AS CPT 


- q+ q- 


0.0 


0.04 ±0.13 


A(7 CPT 


= C £+,K L ~ C £+,K S 


0.0 


-0.04 ±0.15 


AC CPT 


~ ^£+,K L ~ ^£+,K s 


0.0 


-0.06 ±0.13 




q+ 

°£+,K s 


0.7 


0.92 ±0.10 




S £+,K S 


-0.7 


-0.70 ±0.06 




^£+,K s 


0.0 


0.08 ±0.07 




C £+,K S 


0.0 


0.06 ± 0.06 



Table 4. Definition of the T-, CP-, and CPT-asymmetry parameters. These parameters are defined 
as the differences between the S„ o, C^ a coefficients for two reference Ai-exchanged processes and 
those of the corresponding symmetry-transformed transitions. In the central column we show 
the expected values based on the SM CP violation studies at B factory experiments, as given 
in Eqs. (3.3) and (3.4). The right column reports the fit results from one of the 350 simulated 
experiments described in Sec. 4. 



where the a index denotes the other flavor state to that given by a and u a represents the 
fraction of flavor states reconstructed as £~X being a B° (or as £ + X being a B°), with 
3^ (At) given by Eq. (3.2). 

The experiments have employed flavor identification algorithms that analyze tracks not 
associated to the completely reconstructed CP-eigenstate to assign a given category. These 
assignments are determined from different inclusive or semi-inclusive signatures, such as 
isolated primary leptons, kaons, and pions from B decays to final states containing D* 
mesons, and high momentum charged particles from B decays. For this study we use 
different categories with efficiencies e (mistags uj) ranging between 9% and 17% (3% and 
42%), with a total effective flavor identification efficiency of Q = e(l — 2w) 2 ~ 31%. 

The resolution in the reconstruction of the decay time difference by the experiment 



T invariance 


CP invariance 


CPT invariance 


AS+ = 


AS+p = 


ZAO/~iprp U 


AS" = 


A5 CP = 


iA*_)*"iprp U 


A5 CP = A5 CPT 


Z-AOnp — ZAOpi-prp 


Z-XOnn L\iO np 


AS CP = A5 CPT 


ZAOrp Z_iO/~iprp 


ZAOrp — ZAOpip 


AC+ = 


AC+p = 


ZAOz-iprp = U 


AC T =0 


AC CP = 


ZAOpipnp = U 


AC CP = AC CPT 


ZAOrp == ZA(_'/-iprp 


ZAOrp == ZAOpp 


AC CP = AC CPT 


ZAOrp == ZAOz-iprp 


-LAOnp = ZAOpp 



Table 5. Expected values and relations among the asymmetry parameters under invariance of one 
of the three discrete space-time symmetry transformations. 



introduces a smearing such that the observed At might take negative values. Therefore, 
the sign of At, which indicates if the decay to the flavor final state occurred before or after 
the CP final state (as described in Sec. 3), cannot be used directly to disentangle between 
+ and — events. To overcome this problem, we describe the observed decay rate as a 
function of the reconstructed At, At rec , as 



H a ,p(At iec ) ex h+ /3 (At)H(At)®n(St;a Ml J + 

h~ p {-At)H{-At) ® TZ(5t; a Atrcc 



(4.2) 



where H(At), H(—At) are Heaviside step functions, h^JAt), h~a(—At) are given by 
Eq. (4.1), and the symbol <g> indicates the mathematical convolution with the resolution 
function TZ(5t; a At rcc ) , where St = At TCC — At and o"At rcc is the estimate of the At roc un- 
certainty obtained by the reconstruction algorithms. The first term of the sum is related 
to At > (i.e., the tagging B, identified as l + X or £~X, decayed before the other B, 
reconstructed as ccK s or J/ipK L state), while the second term is related to At < (the 
neutral- B meson decayed to a flavor state later). With this construction the distribution 
H. a ^{At rcc ) for At rcc > will be dominated by At > events, but will also contain events 
with At < due to the limited At resolution. Similarly, the distribution T-L a ^{At rec ) for 
At rcc < will contain predominantly events having At < with contribution from At > 
events. The need of distinguishing between At < and At > in the presence of At 
resolution is one of the main complications of this study in comparison to the standard 
CP violation analyses at B factories. This is the reason of the Heaviside step function in 
Eq. (4.2) and the ± index in the upper part of the C a a and S a „ parameters. 

The resolution function 1Z(5t;aAt rcc ) is modeled by the sum of three Gaussian func- 
tions [18] with means and widths proportional to (TAt ICC - The mean scale factor for the 
main component (which includes about 90% of the events) ranges between for flavor 
identification categories containing isolated primary leptons to —0.2 for those categories 



10 



without leptons, to account for the typical offset in the resolution function observed by the 
experiments. 

4.2 Sensitivity results 

We generate a total of 350 experiments containing all 8 pairs of flavor-CP events, each 
of which is then fitted using an unbinned maximum likelihood procedure with the same 
PDF as that employed for the generation of the samples. Each fit determines the best 
8 asymmetry parameters of Table 4, with all other parameters of the PDF, like those 
describing the resolution function, mistags and backgrounds, kept fixed. This is a realistic 
scenario for an experimental analysis since all these parameters are usually determined 
by the experiments using either control samples or detailed simulations of the detectors. 
In order to increase sensitivity and account for the significantly different reconstruction 
efficiencies between the different CP final states (especially between J/ipK L and other 
events) the normalization of the PDF in the fit is performed separately for the different 
CP final states, but simultaneously for £ + X, £~X, At rec > and Ai rcc < events. 

Table 6 summarizes the standard deviation (root mean square, or r.m.s) of the resid- 
ual distribution (fitted minus generated values), the mean of the fit uncertainty, and the 
standard deviation of the fit uncertainty, for each of the asymmetry parameters from the 
350 experiments. The expected sensitivities for the T-violating parameters are 0.15, 0.20, 
0.14 and 0.15, for AS^, AS^, ACp, and AC^, respectively. While the sensitivity for 
AC-p and AC^ is similar, it is significantly worse for AS^, in comparison to ASp", since 
the associated decay rates involve time-dependent distributions with larger tails. The un- 
certainties obtained directly by the fit (parabolic errors) provide a good estimator of the 
true resolution given by the r.m.s of the residual distributions, while no statistically signif- 
icant biases are observed, concluding that all asymmetry parameters have good Gaussian 
behavior. To illustrate the case of an actual experimental analysis, in Table 4 we give 
the complete fit results from one of the simulation experiments. As it can be observed, 
all asymmetry parameters are consistent with the generated values (indicated in the same 
Table), with fit errors in good agreement with the expected sensitivities. 

4.3 Asymmetries 

The difference in the rates for any pair of symmetry-transformed transitions normalized to 
their sum is usually used as observable to probe the symmetry violation. In practice, we 
construct raw asymmetries using the number of events for each pair of transitions in bins 
of Ai + = |At rec |, normalized to the total number of events of the given subsample. This 
normalization is particularly relevant for T (and CPT) asymmetries since it involves com- 
parison of ccK s and J '/ipK L states, which overall have different reconstruction efficiencies. 
Since for a given discrete symmetry there are four possible comparisons between conju- 
gated processes, we have four independent asymmetries. For example, for the B° — > B— 
reference transition the raw T asymmetry would explicitly be defined as 

%- K (Ai+) - H± K (At+) 
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Parameter 


r.m.s residual 


mean error 


r.m.s error 


AS+ 


0.149 


0.148 


0.009 


ASt 


0.201 


0.201 


0.009 


AC+ 


0.139 


0.134 


0.008 


ACt 


0.150 


0.141 


0.009 


AS+p 


0.120 


0.106 


0.005 


A5cp 


0.123 


0.120 


0.005 


AC+ P 


0.103 


0.094 


0.004 


ACcp 


0.103 


0.098 


0.003 


A <?+ 


0.209 


0.213 


0.017 


A6 CPT 


0.137 


0.131 


0.009 


AC CPT 


0.143 


0.147 


0.009 


AC CPT 


0.128 


0.128 


0.007 




0.101 


0.090 


0.004 


5 *+,tf s 


0.060 


0.057 


0.003 




0.069 


0.066 


0.002 


C e+,K S 


0.058 


0.057 


0.002 



Table 6. Root mean square (r.m.s.) of the residual distribution (fit minus generated values), mean 
of the fit uncertainty, and the r.m.s of the fit uncertainty, for each of the asymmetry parameters 
from the 350 simulation experiments. 



where % a JAt + ) = l-L a: p(±At vcc )H(At vcc ). Neglecting mistag and proper-time resolution 
effects, and assuming Eqs. (3.3) and (3.4) in the denominator for simplicity, the raw T- 
asymmetry becomes proportional to (AC™, A/SjT), 



AC A5 
A T (At+) « — 7^cosAmA£ + -\ ^- sin Am At. 



(4.4) 



2 2 

The other three T asymmetries are constructed similarly using the association between 
transitions and time-ordered decay products for reference and T-conjugate processes given 
in Table 1, and are proportional to (AC^, AS^), (AC^ P — ACq PT , AS*q P — AS*q PT ), and 
(ACqp — ACq PT , ASqj, — ASq PT ), respectively. The raw CP and CPT asymmetries are 
constructed following a similar procedure. 

It should be noted, that with the methodology proposed, these asymmetries are only 
used with the purpose of illustrating the symmetry violation effect, through direct com- 
parisons of the time-dependent raw asymmetries from data and the projections of the best 
fit results to the decay rates when we allow for both non-invariance and invariance under 
the symmetry transformation. 
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Figure 2. The four independent raw T- asymmetries corresponding to the four possible compar- 
isons between T-conjugated and reference transitions a) B° — > B- (£ + X,ccK s ), b) B + — > B° 
(ccK s ,£+X), c) B° -> B+ (l + X,J/ipK L ), and d) B_ ->■ B° (J/^K L ,£+X) of Table 1, for one of 
the 350 simulation experiments and combining flavor categories with low mistag (isolated leptons 
and kaons), for a signal enriched region. The points with error bars represent the simulated data, 
the solid (red) curves represent the projections of the best fit results, while the dashed (blue) curves 
represent the projection of the best fit assuming T invariance. 



Figure 2 shows the four independent raw T-asymmetries from data for one of the 350 
simulation experiments, overlaid with the projection of the best fit results together with 
the projection of the best fit under the assumption of T invariance. The T invariance 
hypothesis implies to impose to our most general model in Eq. (3.2) used in the former 
case the 8 restrictions given in the left column of Table 5. This reduces the number of signal 
parameters left free in the fit from 16 to 8. The non-flat shape of the asymmetry curve with 
T invariance imposed is due to the offset of the resolution function discussed in Sec. 4.1. 



4.4 Significance of results 

The significance of T non-invariance can be evaluated by the variation of likelihood (sj.) in 
the space of asymmetry parameters. This variation is determined by the difference between 
the negative log-likelihood of the best fit with the most general model (ln£) and that of the 
best fit imposing the 8 restrictions related to T invariance listed in Table 5 (T-invariance 
point, ItvCt), 



„2 

7 n 



-2Aln£ T = -2(ln£ T -ln£) 



(4.5) 



In the limit that the likelihood function takes a Gaussian shape in the space of asymmetry 
parameters, which is a good approximation in our case as discussed in Sec. 4.2 and sum- 
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marized in Table 6, s 2 -, is properly described by a x 2 distribution with v = 8 degrees of 
freedom. The confidence level (C.L.) can then be obtained by computing the cumulative 
X 2 probability that the value Sj, is exceeded [8] . 

For the same simulation experiment shown in Figure 2 the variation Sj< is 353 units, 
which corresponds to a C.L. equivalent to about 18 standard deviations (here we have 
adopted the convention that 1, 2, 3, ... standard deviations in 8 dimensions have 1 — C.L. = 
0.3173, 4.55 x 1CT 2 , 2.7 x 10" 3 ,...). This is the figure of merit that establishes the degree of 
inconsistency between the best fit result and the T invariance point, i.e., the significance 
that any of the 8 conditions in the left column of Table 5 is not satisfied, or equivalently, 
how significant is the combined difference between the solid (red) and dashed (blue) curves 
for the four asymmetries in Figure 2. 

The significance of CP and CPT non-invariance can be evaluated similarly. For the 
same experiment as previously, the variations in s 2 for the CP- and CPT-invariance points 
are found 536 and 2 units, corresponding to about 22 and less than 1 standard deviations, 
respectively. These values are consistent with T non-invariance compensating CP viola- 
tion so that CPT remains invariant, in agreement with the input parameters used in the 
generation of the simulated experiment. 

5 Conclusions 

This work concerns the study of microscopic Time Reversal Violation in the fundamental 
laws of physics. The observed time asymmetries in current macroscopic and microscopic 
phenomena are not connected to this problem and they can occur in theories with exact 
T invariance. 

A direct evidence for T violation means an experiment that, considered by itself, 
clearly demonstrates T violation independent of the results for CP violation. For transition 
processes, the antiunitarity of the operator implementing the symmetry transformation 
implies that a genuine test of T violation needs an asymmetry under the interchange of in 
and out states in the dynamical evolution of the system. Nowadays, there is no experimental 
result providing direct evidence of genuine T violation. The measured asymmetry among 
the probabilities K° — > K° and K° — > K° cannot be interpreted as such since, being 
CPT-even transitions, CP and T are experimentally identical, no matter whether there 
is CPT invariance or not. Other experiments that could provide evidence involve a non- 
vanishing expected value of T-odd observables for stationary, non-degenerate states, like 
electric dipole moments, not yet observed. For unstable particles, where we might expect 
large effects, the associated irreversibility seems to prevent a direct test of T symmetry 
through the exchange of in and out states. 

This paper discusses how to overcome this irreversibility problem and defines the pre- 
cise steps for implementing an actual experiment able to obtain direct evidence of genuine 
T violation in the time evolution of a neutral- B meson, taking place between the two decays 
of the B°B° system produced at B factories. The essential ingredients are the quantum 
mechanical entanglement imposed by the EPR correlation, in particular, the first decay 
of one B prepares the quantum state of its living partner, and the experimental ability 
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to unfold the ordering and difference of the two decay times. Identifying the B decays 
into definite flavor (£ + X or £~X) and CP-eigenstates {J/ipK s or J/tpK L ), the difference 
between the decay rates for (£ + X, J/tpK s ) and (J/ipK L , £~X) transitions at the same 
decay time difference is our signature of T violation. There are three other independent 
signatures which can be constructed from the eight time-ordered pairs of B decays. These 
proofs are independent of the results for CP violation, as shown by the fact that, given a 
reference (£ + X, J/ipK s ) for example, the T-, CP-, CPT-transformed, and At-exchanged 
processes are all experimentally different. 

To demonstrate the feasibility of the proposed methodology and how it can be applied 
for a real data experimental analysis, we have performed a realistic Monte Carlo simulation 
of data samples with similar size and properties to those currently available at B factories, 
taking into account misidentification of CP eigenstates, mistakes in the flavor identification 
and resolution in the reconstruction of the decay time difference. The need of resolving 
simultaneously the ordering and difference of the two decay times in the presence of reso- 
lution is one of the main challenges in comparison to standard CP violation studies. For 
the purpose of illustrating the T-violating effect, we have presented a direct comparison of 
the time-dependent raw asymmetries from data and the projections of the best fit results 
to the decay rates with and without T non-invariance. 

Using a large number of generated experiments, we have inferred the expected sensi- 
tivities of the asymmetries in the parameters which determine the time-dependent decay 
rates, together with the expected significance of T non-invariance. A large significance, 
equivalent to about 18 standard deviations, is foreseen with the currently available data, 
providing the expectation of a discovery of direct, genuine Time Reversal Violation in the 
time evolution of a neutral- B meson. 
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A The |5_) and \B+) states 

The I-B-), \B+) states are experimentally identified as those filtered by the observation 
of the decay to definite CP eigenstates. The only requirement needed for the analysis 
is their existence, independent of the underlying theory (CPT invariant or not), with the 
ortogonality property {B_\B + ) = 0. In this Appendix we construct these states explicitely. 
Neutral kaons decaying inside the geometrical acceptance of detectors surrounding 
their production point are usually reconstructed through their decay to two pions. Thus 
we can first define \BJ) as the state filtered by the decay into J/ipK + , K + being the 
neutral--£f state decaying to tttt, a pure CP-odd state. The state \B + ) is then defined as 



15 



the state orthogonal to |-B_), (B+\B—) = 0, which cannot decay into the J/ipmr final state, 
{J/i/j7m\T\B + ) = 0, where T is the transition operator. The state \B + ) can be written in 
terms of flavor eigenstates as 

\B + )^^ + [\B°)-a\B°)], (A.l) 

where a = = — and N_i_ is a normalization constant. Since \B_) is orthogonal 

(J/^Tnr\T\B ) + 1/5 

to |-B+) it then follows 

\BJ) =N_ [\B°) + a*\B )] . (A.2) 

Analogously, we define \B + ) as the state filtered by the decay into J/ipK L , a CP-even 
state up to 0(1O~ 3 ) due to CP violation in the neutral- kaon system. We note that we could 
have defined \B + ) through its decay to J/ipK_, K_ — > 7r 7r 7r , in this case a pure CP-even 
state, but this final state cannot be reconstructed since long-lived neutral kaons tend to 
interact hadronically inside the detectors before they can undergo decay. Of course, both 
definitions become operationally identical when CP violation in neutral-kaons is neglected. 
The state \B—), defined as its orthogonal state, (B-\B+) = and (J /ipK L \T\B-) = 0, is 

|,B_> =N_ [|i?°> -/3|^°>] , (A.3) 

(J/ibK \T\B°) 

where B = = — . Therefore the state \B+) can be explicitly written as 

(J/if>K L \T\B°) 

\B + )=N + [\B°)+P*\B )}. (A.4) 

Let us note that here we keep separate the definitions of the states |-B-) and \B + ), 
which are observed through their decays to the J/ipinr and J/ipK L final states, from the 
states I-B-) and \B + ), produced exploiting the EPR correlations in the entangled BB 
system. The two bases are the same when the orthogonality condition (B_\B + ) = is 
fulfilled. As a and (3 correspond to opposite CP final states, we have a/3* = — 1 when only 
one weak amplitude is responsible of the decay (\a\ = \(3\ = 1). Consequently as all these 
considerations apply to our experimental framework, 

(-B-I-B+) = 1 + a/T = , (A.5) 

which was to be demonstrated. 
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